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ABSTRACT 

We present theoretical predictions for photometric and spectroscopic signa- 
tures of rings around transiting extrasolar planets. On the basis of a general 
formulation for the transiting signature in the stellar light curve and the velocity 
anomaly due to the Rossiter effect, we compute the expected signals analytically 
for a face-on ring system, and numerically for more general configurations. We 
study the detectability of a ring around a transiting planet located at a = 3AU 
for a variety of obliquity and azimuthal angles, and find that it is possible to de- 
tect the ring signature both photometrically and spectroscopically unless the ring 
is almost edge-on (i.e., the obliquity angle of the ring 9 is much less than unity). 
We also consider the detectability of planetary rings around a close-in planet, 
HD 209458b {9 ^ 90° - iorbit ~ 3°. 32), and Saturn {9 ^ 26°. 7) as illustrative 
examples. While the former is difficult to detect with the current precision (pho- 
tometric precision of 10~^ and radial velocity precision of 1 m s~^), a marginal 
detection of the latter is possible photometrically. If the future precision of the 
radial velocity measurement reaches even below 0.1 m s~^, they will be even 
detectable from the ground-based spectroscopic observations. 
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Introduction 



In 1903, Hantaro Nagaoka, a professor in pliysics, Imperial University of Tokyo, proposed 
an atomic model which consists of a large number of particles of equal ma ss arranged in a 



circle at equal angular intervals around a central particle of large mass (INagaoka I Il903 
19041 ) ■ He referred to this model as a Saturnian mo del, and published the hypothesis prior 



to the famous Rutherford model (IRutherford Ill91ll ). often referred to as a planetary model 
of the atom. This remarkable example in history illustrates how our Solar system, and even a 
Saturnian ring, inspired the best physicists in the world to contemplate the quantum world. 

Thus it is interesting to contemplate how the properties of more than 300 extrasolar 
planets discovered so far beyond our Solar system will inspire future scientists. Among 
them, 52 systems are known to exhibit transiting signatures. Detailed photometric and 
spectroscopic observations of such transiting planets provide unique opportunities to look 
into the nature of extrasolar planets, which is otherwise unaccessible including the planetary 
radius, average density, and atmospheric composition. Ongoing/upcoming missions such as 
COROT and Kepler will significantly increase the number of transiting planets, and even a 
wider range of research will be made possible in the next several years. One specific example 
i s the exploration of the spin-orbit rnisalignment of extrasolar planet using the Rossite r effect 
feossiter lE^ : lMcLaughlin1ll924j : IPetrie Ill938] : iKopal Ill94i Il945l : Iflosokawa lEgSsh : when 
a planet transits its host star, the light coming from the star is partially blocked. Photomet- 
rically this produces a characteristic dimming of the stellar light curve, and spectroscopically, 
an asymmetric distortion of the absorption line profiles of the stellar spectrum due to the 
occultation of some of the velocity components of the stellar absorption lin es. Thus th e 



latter yields an apparent anomaly in the radial velocity curve of the host star (iQueloz et al. 
2000h . 



In a previous paper (lOhta. Taruya. fc Sutdl2005l . hereafter Paper I) , we derived analytic 
templates for the v elocity anomaly using a perturbative expansion. Motivated by that work, 
Winn et al.l (120051 ) combined the best observational datasets currently available for the HD 
209458 system, and detected for the first time a small misalignment (A = —4°. 4 ± 1°.4) 
between th e spin axis of the central star and the orbital axis of the planet. Subsequently 
Wolf et al.l (120061 ) applied the analytic formulae of Paper I, and obtairied th e limit, A = 
+ ir ± 14°, for the anomalously dense planet, HD 149026 JSato et al.ll2005h . Also, the 



spin-orbit alignment of the HD 189733 system were me asured bv iWinn et al. 
very small misalignment of A = — 1°.4 ± l°.l was found (IWinn et al. 



2007, 



aa 



J2006h and 



Narita et al. 



20071 . see also). Recently measurements of spin-orbit alignment have been successful for a 
number of transiting planets. It is just a matter of time to obtain the statistical distribution 
of the misalignment angle A, which provides a unique probe of the origin and evolution of 



- 3 - 



the angular momentum of the planetary systems. 



In the present paper, we consider an even more ambitious possibility to detect rings 
around transiting planets. It is interesting to note that the rings aro und Uranus and Neptune 
were i ndeed discovered duri n g their occultation of b ackground stars (lElliot. Dunham fc Mink 



1977; Hubbard et al. 1977 



19861 : lEsposito 1 12006| ) . While it is not clear whether or not 



extrasolar planets have rings and/or satellites as those in our solar system, their detection, 
if successful, would definitely mark an important milestone in astronomy and planetary 
sciences; transiting planets are the unique targets for that purpose. From the methodology 
developed in this paper, one can further determine the orientation of the spin of the orbiting 
planet under a reasonable assumption that the orientation of the ring is aligned to the 
planetary spin axis; this is approximately the case for planetary rings in our Solar system. 



Indeed we are not the first to explore the detectability in detail. iBrown et al.l (120011 ) . for 
instance, put photometric constraints on possible satellites and rings around the first transit- 
ing planet, HD 209458b. They assume that an optically thick two-dimensional ring extending 
from the surface of the planet up to the outer radius -Rout- Hot Jupiters like HD 209458b 
are supposed to be roughly tidally locked to the host star. Thus they also assume that the 
ring plane coincides with the planetary orbital plane, and found that i?out should be less 
than 1.8 times the planetary radius. Given the approximate edge- on nature of the as s umed 
ring, the upper limit is fairly strong and interesting. More recently iBarnes fc Fortneyl (l2004j . 
hereafter BF04) presented a light curve of ringed planets over a wide range of the parameter 
space, and discussed the photometric detectability of the additional ring signature. 

Our present paper differs from the previous work in that we examine the spectroscopic 
detectability of planetary rings using the Rossiter effect in addition to the photometric one. 
Given the current observational sensitivities, the photometric signatures are easier to detect 
than the spectroscopic ones in most cases. Nevertheless, the complementary nature of the 
spectroscopic detection/confirmation of rings is important because the signature of rings is 
very subtle and may be mimicked/disguised by unknown noise and systematics. Moreover 
the spectroscopic signature may be more pronounced for transiting planets in rapidly rotating 
stars because the velocity anomaly due to the Rossiter effect is simply proportional to the 
stellar spin velocity. 

We also present improved and more accurate analytic formulae of the photometric light 
curve and the spectroscopic velocity anomaly for a transiting planet without a ring than 
those in Paper I. While those for a planet with a ring need to be computed numerically in 
general, such analytic formulae are very useful in performing the parameter fit for a ringless 
planet model, and also in checking the accuracy of the numerical results for a planet with a 
ring. 
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The plan of the paper is as follows; we describe our basic assumptions for the planet and 
ring system in §2. Section 3 summarizes the general formulation to compute the photometric 
and spectroscopic signals for the planet and ring system. Analytical results for the face-on 
ring are presented in §4, which are in fact the improved version of the results of Paper 
I. Numerical results for an arbitrarily oriented ring are shown in §5, and we discuss its 
detectability in detail. Finally §6 is devoted to conclusions and imphcations of the paper. 



Let us begin with describing several definitions and assumptions concerning the geomet- 
rical configuration of the planetary system and the properties of the host star, the orbiting 
planet and a surrounding ring. Table [1] lists the notations and the parameters of our model, 
which are used in the analysis below. 



Fig. 1. — Planetary transit accompanied by a ring. The internal parameters of the star 
and the planet are shown: stellar radius Rs] projected angular velocity Qs sinJ^; impact 
parameter of planetary orbit, b; angle between the x-axis and the orbital plane. A, which 
represents the spin-orbit alignment projected onto the x — z plane. 



We assume that the star (radius Rg) and the planet (radius Rp) follow the exact two- 
body Kepler orbit, ignoring perturbations from possible outer planets. We adopt a coordinate 



2. Model assumptions 



z 




2.1. Star-planet system 
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system centered at the star so that its y-axis is directed toward the observer, and the z-axis 
is chosen so that the stellar rotation axis lies on the y-z plane. We assume rigid-body stellar 
rotation with VLg being its spin angular velocity. In this coordinate system, the components 
of the stellar spin angular velocity reduces to 

f2s = (O,fiscos/s,ilssin JJ, (1) 

where Ig is the inclination angle between the stellar spin axis and observer's direction. 

Since we are mainly interested in the transit phase whose duration is much shorter than 
the entire orbital period, the actual eccentric orbit is close to a circular one in the vicinity of 
the planetary transit. Thus we adopt the circular orbit in the present analysis while most of 
our formulae below are applicable to a non-circular orbit by properly mapping the projected 
position of the planet into the observed time. Under the circular orbit approximation, the 
planetary orbit is simply characterized by the impact parameter h and the projected angle A 
between the stellar spin and the planetary orbital axes, in addition to the semi-major axis a 
(Fig. d]). The impact parameter h is equal to 6 = acosiorbit/-Rs, where iorut is the inchnation 
of the orbit. 

We adopt the linear limb darkening alone for the surface intensity of a stellar disk in 
Paper I. In this paper, we take account of the effect up to quadratic order: 

htM = ^0 [1 - - /i) - U2{i - ^lf] . (2) 

The two constants, ui and M2, characterize the amplitude of the limb darkening, and /i is the 
directional cosine between the line of sight and the normal vector to the local stellar surface: 

/, + 

2.2. Planet and ring 

We consider a simple model of ring around a planet as shown in Figure [2l The ring 
is two-dimensional (geometrically thin) and circular. Its inner and outer radii are denoted 
by i?in and i?out which satisfy Rp < i?in < -Rout- We assume that the ring moves along 
the planetary orbit with constant obliquity angles {9, cj)). The light transmitted through the 
ring diminishes due to the absorption and/or the extinction by particles in the ring. We 
introduce the normal optical depth of ring, r, so that the transmitted flux through the ring 
is proportional to exp(— r//5), where (3 is the cosine of the angle between the normal vector 
to the ring surface and y-axis, and given hy (3 = sin 9 cos (f). We adopt r = 1 for simplicity 



(3) 
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(see Table HI); the optical depth of the B ring of Saturn is in the range of r 
jMurrav fc Dermottl1l999l ) . 



0.4 ~ 2.5 



X 




Rout 



Fig. 2. — Schematic configuration of the planet accompanied by a ring. The parameters of 
the planet and the ring include planetary radius Rp, radius of inner edge of ring R-^^, radius 
of outer edge of ring -Rout, obliquity angle 6 (the angle between the normal vector to the ring 
surface and 2;-axis), rotation angle (the angle between the |/-axis and the projected normal 
vector of ring surface onto the y — z plane). 



Barnes fc Fortneyl (120041 ) examined the influence of diffraction on the photometric light 



curves, and found that it may leave a detectable signature depending on the size of ring 
particles. We ignore the diffraction of the ring particles here for two reasons; first, our 
primary purpose is to develop a methodology which directly tests the presence of rings, 
and second, the diffraction effect itself does not exceed the primary effect of the star-light 
dimming. In any case the effect of diffraction can be easily incorporated into our numerical 
scheme if the ring properties are specified. 



3. Photometric and spectroscopic transit signatures for a planet with a ring 

3.1. Photometric signals 

The photometric light c u rve fo r transiting extrasolar planet with a ring has been com- 
puted by iBarnes fc Fortneyl (120041 ) in detail. For definiteness and later convenience, we 
summarize their basic equations in terms of our own notations here. 



The dimming of the star light during the transit phase is quantified by evaluating the 
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fraction of star light blocked by the planet and/or intercepted the ring. We define the relative 
fiux ratio F: 



Is{x, z)dxdz 



(4) 



Here, the functions /(x, z) and /s(x, z) represent an intensity at a given position (x, z) of the 
stellar surface. The stellar intensity model ([2]) and the ring model described in §2.21 yield 
the explicit form of the function /(x, z): 



I{x, z) 






-^star(/^) 



region outside the star, i.e., + z"^ > 
region inside the planet, i.e., (x — Xp)'^ + {z ^ ^•.p 

region transmitted through the ring, see conditions (j6]) ' 
otherwise 



where the direction cosine /i is given by equation ([3]). The region where the light coming 
from the star is transmitted through the ring satisfies the three conditions: 



R^r. < 



Zp) sin I 



COS( 



+ 



sin 9 



^ d2 
^ -"-out' 



z'<Rt 



X 

(x-Xp)' + iz~Zpf>Rl 



(6) 



The quantities, Xp and Zp, represent the position centered at the planet, which manifestly 
depend on time. 



The stellar intensity outside the transit phase, Is{x,z), is simply written as 

^""^^'^^ { /star(/i) : msiae X 

, the denominator of equation 



outside the star i.e., x^ + > R^ 
inside the star i.e., x^ + z"^ < R^ 



Substituting the stellar intensity model ([ 
grated analytically: 

j j Is{x, z)dxdz = t^IqRI ~ y ^ y 
Note that this is independent of time. 



(7) 

can be inte- 

(8) 



3.2. Spectroscopic signals 



While planetary rings may be detectable with photometric data analysis alone, the 
spectroscopic data provide an independent and complementary check of the presence of 
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rings. Since the two methods would suffer from very different noises and systematics, the 
agreement between them significantly increases the credibility of the positive detection. 



The spectroscopic detection method relies on the Rossiter effect ( iRossiter Ill924l ) which 
represents the "apparent" radial velocity modulation arising from the spectral line distortion 
due to an occultation of a part of the rotating star. During the passage of the transiting 
planet, the light from the stellar surface is asymmetrically blocked off. As a result, the 
line-profile-weighted mean position for a specific absorption or emission line taking account 
of the stellar rotation is apparently shifted. In Paper I, we derived analytical expressions for 
radial velocity curves for extrasolar transiting planets for the first time, but the effect of rings 
was ignored. The quantitative modelling of the photometric and spectroscopic signatures of 
rings will be considered in detail below. 

Formally the radial velocity anomaly due to the Rossiter effect is expressed as 



X /(x, z)dxdz 

Avs = -^ssinls"^^^ . (9) 

I{x, z)dxdz 



The above expression, which relates the velocity anomaly to the first moment of the line 
profile, is applicable also to transiting planets with a ring if equation ([5]) is substituted 
into I{x,z). In general, the above integral cannot be performed analytically but needs 
to be done numer ically (except for a face-on ring as described in §1]). We note here that 



Winn et al.l (120051 ) pointed out the presence of ~ 10 percent systematic difference between 
the perturbative results of Paper I and their numerical fits to the outputs of the data analysis 
pipeline. This is supposed to come from the fact that the approximation of the first moment 
of the line profile as Avs is not sufficiently accurate. The systematic difference should be 
sensitive to the specific data analysis routine, and hard to evaluate generally. Moreover it 
is unlikely to affect our conclusions which depend on the difference of Avg of planets with 
and without rings. This is why we use equation ([9]) to evaluate the Rossiter effect in what 
follows. 

In practice, we proceed as follows; we first calculate the integral over an entire surface of 
the stellar disk, which can be done analytically. Then we numerically evaluate the integrals 
over the regions overlapped with the planet and the ring, and subtract these from the analytic 
result without a planet nor a ring. In order to cope with the complicated geometry of the 
integration, we set square grids (256 x 256 cells) enclosing both the planet and the ring with 
a length of 2i?out on a side. Judging from the location of the center of each cell relative to the 
stellar surface, we assign the intensity /(x, z) in equation ([5]), and compute the contribution 
from the cell. We test the convergence and the accuracy of our numerical integrations by 
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comparing with the results with finer grids (512 x 512 cells). We find that the fractional 
difference of the photometric signal is less than 3 x 10~^, and that the spectroscopic signals 
agree within 0.08ms~^, both of which are negligibly small for our current purpose. 



4. Analytic formulae for a planet with a face-on ring 



Equations (jlj) and iQ can be integrated (almost) analytically for a planet with a face-on 
ring, i.e., {6, (p) = (90°, 0°). These analytic formulae help our understanding of the basic 
observational signatures of a ring from the photometric and the spectroscopic data. They 
are also useful in making sure of the accuracy of the numerical results for more general cases. 

We find that equation (jl]) reduces to 

^ _ 7r(l - ui/3 - U2/6) - A{p, 7ps) - (1 - e-") {A{p, 70,) - A{p, 7;,)} 

7r(l -M1/3-U2/6) ^ ' 

for the relative fiux ratio, and equation iQ to 

5(p,7ps) + (1 - e-") {5(p,7os) - S(p,7i.)} 



Avs = Xr, Vts sin Is — , , , , , , , X r . / 

' 77(1-^/1/3 -n2/6)-A(p,7p,)- (l-e-^){A(p,7os)- A(p,7i3)} 



:ii) 



for the radial velocity shift, where 7ps = Rp/Rg, 7is = Rin/Rs, 7os = Rout/Rs, and p = 
^yx^ + Z^/Rg. The quantities, A and B, are given in terms of the normalized planet 
position p as follows: 



Mp, 7) 



f 7772 [l-Ui- U2 (2 - p2 _ y/2) + {Ui + 2U2)Wi] 

(1 - ui - 3^2/2) [72 cos"^ (C/7) + sin"^ zq - pzo] 
+ {u2/2) p [(p + 2C)7' cos-i (C/7) - (pC + 27')] 

+ {UI + 2U2)W3 



p < 1 -7 



1 - 7 < P < 1 + 7 (12) 



1 +7 < P 



and 



' 717^ [1 - Ml - U2 (2 - p2 - 7^) + (mi + 2u2)W2] : p < 1 - 7 



B{p, 7) = < 



U2 (2 - p2 - 7^)] [7^ cos-i (C/7) - ^oC] 

-{Au2/i)pzl + {Ui + 2M2)W^4/P 



: 1 — 7<p<l + 7 



1 + 7 < P 



(13) 
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The variables zq and C are also written in terms of 7 and r/p = p — 1 as 



^0(7) 



(7^ - V') iVp + - 7^ 



C(7) 



Vl + + 7^ 



(14) 



The functions VFj (i = 1 ~ 4) in the above expressions are in fact analytically intractable, 
but we find accurate analytical fitting formulae which are explicitly given in Appendix |Al 

If we set r = 0, the above analytic expressions reduce to a case of a planet with- 
out a ring. However they a re superior to, and thus should replace, our previous results 
(lOhta. Taruya. &: Sutdl2005l ) in two ways; first, we now consider a quadratic limb darkening 
(eq.[2]), instead of a linear limb darkening {u2 = 0), and second, the approximate expressions 
for Wi given in the appendix are more accurate by a factor of a few than the previous one; for 
typical parameters 7ps < 7is < 7os ~ 0.2, the accuracy is typically within a few percent level. 
We use the above formulae extensively when we discuss the effect of rings in the following 
analysis. 



5. Predicted photometric and spectroscopic signals 

Even in the simple model that we assume here, the photometric light curve and the 
spectroscopic velocity curve are characterized by 12 and 14 parameters, respectively. The 
former includes the radii of the star and the planet {Rg and Rp), the inner and the outer 
edges of the ring (_Rin and -Rout), the limb-darkening parameters {ui and ^2), the orbital 
parameters of the system (semi- major axis a, inclination of the orbital plane iorbit, and the 
orbital rotation period Porbit), the obliquity angles of the ring (6* and 0), and finally the 
normal optical depth of the ring, r. The additional two parameters for the latter are the 
projected stellar surface velocity Vsinlg and the projected misalignment angle A between 
the stellar spin and the orbital plane axes. 

In order to avoid the tidal disruption, the outer radius of the ring should not exceed the 
Hill radius: 

< . f^V".. 15.«. f#V" f#V" (^) . (15) 



3MsJ \MjJ \MsJ V0.05AU 



Gaudi. Chang fc Han I (120031 ) discussed further theoretical constraints on planetary rings. 



Planetary rings consisting of ice sublimation if the planetary semi-major axis satisfies 

r \ / r \ 1/2 
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where Tgub = 170K is the subhmation temperature of water ice, and is the luminosity of 
the host star. For rocky rings, the Poynting Robertson drag force, viscous friction due to the 
planetary exosphere, torques due to (shepherd) satellites, and internal scattering become 
important. For close-in planets, the Poynting Robertson drag time-scale tpR is generally 
shorter than the viscou s friction decay time-scale un less the density of the exosphere is 



large, pes ~ 10 ^^gcm ^ (IGoldreich fc Tremaine Ill982l ). The former is given as 



p r / n \'i 



tpR-lOV^^^ hrTTTT , (17) 
Igcm 1cm VO.IAU/ 

where p and r denote the density and the radius of ring particles. Thus rings around close- 
in planets would not live longer than t-pK unless any stabilization mechanism like shepherd 
satellites operates. Thus it is not clear if close-in planets are accompanied by rings. 

Based on these estimates, we adopt a fiducial set of parameters listed in Table [T], which 
is motivated by the existing close-in planetary system HD 209458, but with different semi- 
major axis, a = 3AU. As illustrative examples, we also consider the close-in planet and the 
Saturn system as a member of solar-system planet. We first present a special case of the face- 
on ring (that can be analytically computed) so as to elucidate the basic features ( §5.11) . and 
to quantify the parameter dependence of the photometric and spectroscopic signals in §5.21 
Next in §5.31 we discuss how to identify the signature of rings by looking for the departure 
from the fit to a planet model without ring. Based on these discussions, the detectability of 
a ring system is addressed in §5.41 in addition to future prospects for constraining the ring 
parameters. 



5.1. Signature of the face-on ring 

Figure [3] shows the velocity anomaly by the Rossiter effect (upper) and the light curve 
of photometry (lower) for a planet with a face-on ring, i.e., 6 = ip = 0. The central transit 
time is chosen as t = 0. These curves are plotted using the analytical expressions ffTOl) and 

(ED. 

Compared to the curves of a planet without a ring (dashed lines), the planetary ring 
produces larger radial velocity shifts and flux dimming (solid lines). Also the transit duration 
increases. As we will show below, however, most of those features due to the ring resemble 
that of a ringless planet with a larger radius. The vertical dotted lines in Figure [3] indicate the 
epochs when the boundary edges of the planet and the ring contact the edge of the stellar 
disk (i.e., the ingress and egress phases). These somewhat discontinuous features in the 
derivatives of the curves, if detectable directly at all, would be an unambiguous conflrmation 
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Fig. 3. — Schematic illustration for photometric and spectroscopic signatures of transiting 
planet with and without a ring (see Tabled] for the fiducial values of the model parameters). 
The top panel shows the velocity anomaly curve for the Rossiter effect. The ring is assumed 
to be face-on, i.e., 6 = 90°, = 0°. The curve is analytically calculated according to the 
expression (fTTl) . The bottom panel shows the light curve during transit. The ring is also 
face-on. The curve is analytically calculated according to the expression (|TOi) . 



of the presence of the ring. In reality, however, it may be difficult to locate the features from 
the sampled data with a finite exposure time. 

A more realistic methodology is to first search for the best-fit parameters assuming no 
ring at all, and then to loo k for any systematic rin g signatures in the residuals between the 
observed data and the fits (IBarnes fc Fortneyl 120041 ). Even so, the detection of rings requires 
better resolutions in sampling time, as well as in photometric and spectroscopic accuracy, 
than those of current data. 



5.2. Dependence on the ring parameters 



In general, the photometric and the spectroscopic signals due to rings need to be eval- 
uated numerically. Here we discuss their sensitivity on the four ring parameters {6, (p, Rm 
and -Rout)- Figure m shows some examples with slightly different parameters from our fidu- 
cial values; ^ = 45° ^ 70° {dotted), = 0°^ 45° {short- dashed), Ri^ = l.5Rp -> l.ORp 
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Fig. 4. — Same as Figure [3], but for a ring with different set of parameters. The top panel 
shows the velocity anomaly from Rossiter effect, while bottom panel plots the light curve 
during transit. In each panel, solid lines represent the results for the fiducial model (Table 
[1]). Other lines show the results when one of the fiducial parameters is varied; 6 = 70° 
(dotted), = 45° (short-dashed), i?in = Rp {long- dashed), and -Rout = ^.bRp (dot-dashed). 



(long-dashed) and -Rout = 2.0-Rp 2.5-Rp (dot-dashed). 

In order to clarify how the predicted Avs and F are sensitive to the underlying param- 
eters of the system (9, 0, -Rin and -Rout), we compute their dependence as we did in Figure 
9 of Paper I. Specifically, Figure [5] illustrates the variation of the radial velocity shifts Avs 
with respect to p = -Rin, -Rout and 9: 

5Avs Avs(p + dp;t) - Avs(p;t) 

— = lim — . (18) 

omp dp^o ap/p 

Similarly we define the variation of the photometric light curves 6F: 

^ F(p+dp.t)-F(p;t)_ j^g^ 



S\np dp^o dp/p 

which are also plotted in Figure [5l Since the dependence on 9 is small, we multiply the 
values of 6Avs/6\n9 and 6F/6\n9 by a factor of 10 in the plot. 

Figure H] indicates that a different set of the ring parameters mainly changes the ampli- 
tude of signals during the full transiting phase, and that the bumpy structures around the 
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Fig. 5. — Variation of the photometric (upper) and spectroscopic (lower) signals with respect 
to the ring parameters, i?in, -Rout, and 9. The curves for the latter case are multiplied by a 
factor of 10 for clarity. 



ingress and the egress phases, which indeed are the important signatures of possible rings, 
are not affected much. This can be understood from the weak dependence on the angle 
parameters in Figure [51 If we changes the value of one of the ring parameters from p to p', 
the photometric and spectroscopic signals become 

« (20) 

Ai..(p';i) « Av.(p;t)+(^)^-^. (21) 

\ p J mp 

Those approximate perturbation formulae combined with Figure [5] explain the behavior of 
Figure m fairly well, in particular the weak dependence on 6. 

We note that the three ring parameters {9, R^^ and -Rout) symmetrically change the 
photometric light curves around the central transit time and anti-symmetrically change the 
radial velocity shifts (Fig IS]). Such behavior is quite similar to the variation with respect to 
the planetary and the stellar internal parameters (see Fig.9 of Paper I). Thus, the overall 
change of the amplitude is not a sensitive measure of the presence of a ring. Therefore the 
precise measurement of the bumps around the ingress and the egress phases, even if difficult, 
provides the key in the detection and the characterization of the planetary ring. 
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5.3. Signatures of rings 

The results of the previous subsections now enable us to discuss in detail the signatures 
of a ring system in photometric and spectroscopic data. Given real data set, first thing that 
we should do is to fit the photometric light curve and/or the spectroscopic velocity anomaly 
of a transit planet using a model without ring. If the residuals between the data and the 
best-fit model are significantly larger than the observational errors, one suspects the presence 
of a ring. 

More specifically, let us introduce the above residuals as 

Res[vs] = Avf' - Av^'''-^' (without ring) (22) 

and 

Res[F] = F°^' - F^"""^-^^ (without ring). (23) 



Lovis et al.l (120061 ) has already achieved ~ Im s ^ sensitivity in the radial velocity mea- 



surem ent . Also the HST photometric accuracy is shown to be better than 0.1% (IBrown et al. 



200ll ). which is expected to be the case for upcoming space missions as well. Figure E] illus- 
trates examples for Res[fs] and Res[F]; the upper and lower-left panels show our fiducial ring 
model (Tabled with three different orientations, {9 = 45°, = 0°), {6 = 90°, = 0°) and 
{9 = 90° — Zorbit = 3°. 32, = 0°). The last case corresponds to a close-in planet, HD 209458b, 
with ring. The semi-major axis and inclination angle of the orbital plane with respect to our 
line-of-sight are set to a = 0.0468AU and Zorbit = 86°. 68, respectively. Since close-in planet 
is likely to be a tidally-locked system, we assume that the ring plane lies in the planetary 
orbital plane. Thus, 9 is equal to 90° — iorbit- The lower-right panel is intended to illustrate 
the expected signals for the Saturn system (B-ring) just for comparison; 9 = 26°. 7,0 = 0°, 
i?in = 1.53i?saturn, -Rout = l-95i?satum, ^ siu = l,880ms~^. Note that in the lower panels, 
the residuals obtained from the numerical fit are also plotted as dashed lines, in addition to 
the results from the analytical fit (solid). The residuals, Resff^] and Res[F], clearly indicate 
the characteristic pattern expected for a ring. More importantly, the above signatures are 
indeed at a few a level of the currently best observational sensitivities, both photometrically 
and spectroscopically. Even though it is not clear if those values that we adopted for the 
ring parameters are realistic, the result is certainly encouraging. 

Figure [7] displays the results when one parameter of the ring system is changed from 
the fiducial values, 9 = 45°, = 0°, -Rin = l.5Rp, and i?out = 2.0-Rp (the upper-left panel 
of FigM); 9 = 70° (upper-left), = 45° and Zorbit = 89° .95{upper-nght) , and R^^ = l.ORp 
(lower-left) and -Rout = 2.5-Rp (lower-right). The behavior is fairly easy to understand in 
an intuitive manner. As long as the clear gap between the planetary surface and the inner 
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Fig. 6. — Spectroscopic and photometric signatures of a ring. Each panel displays the velocity 
anomaly due to the Rossiter effect Avg, the spectroscopic residual Res[fs], photometric light 
curve F, and the photometric residual Res[F] from top to bottom. Solid and dashed curves 
in Avs and F present our predictions and the best-fit without ring model. Note that in 
lower panels, the residuals obtained from the numerical fit are also plotted as dashed lines, 
in addition to the one obtained from the analytical fit (solid) ; 6 = 45°, = 0° (upper-left), 
9 = 90°, = 0° (upper-right), tidally-locked close-in planet adopting the orbital parameters 
of HD 209458 system, i.e., a = 0.0468AU, 6 = 3°.32, = 0° (lower-left), and Saturn system 
(B-ring) with e = 26°.7,0 = 0°, i^m = 1.53i?saturn, Rout = 1.95i?saturn, Vsinh = l,880ms-^ 
(lower-right) . 
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Fig. 7. — Same as the upper-left panel of Figure [6l but with a different parameter set; 
6 = 70° (upper-left), = 45° and iorbit = 89°. 95 (upper-right), and Rin = IRp (lower-left) 
and i?out = ^.bRp (lower-right) . 
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edge of the ring exists, a ring produces a larger residual as 9 approaches 90° (face-on); an 
edge-on system is hard to identify. A non-zero value of (j) generates an asymmetry in the 
residuals depending on the value of the inclination angle of the orbit iorbit; Res[F] is not 
symmetric with respect to t ^ —t, and Res[fs] loses a symmetry with respect to the origin. 
The latter effect partially resembles the behavior of the spin-orbit misalignment and results 
in a systematic offset of the value of A. While we restrict our analysis here for A = 0° case, 
this needs to be kept in mind in future precise analyses of ring systems in general. If the 
width of the ring, AR = i?out — Rm, is larger, the signals themselves always increase, but 
the residuals may decrease in some cases (lower-left panel of FigjTj). A trivial example is a 
face-on ring system with Ri^ = Rp which cannot be distinguished from a single planet if the 
ring is opaque. The gap between the planet and the ring leaves an important observational 
signature of the ring. 



Incidentally IColwell et al.l (120071 ) showed that the line-of-sight optical depth of Saturn's 
ring is roughly independent of the viewing angle. If this is the case, we should set (3 = 1 in 
equation ([5]). Thus we repeated the same calculations for a variety of different values of 6 
by setting (3 = 1 instead. We found that the amplitudes of the observable ring signatures, 
i.e., residuals defined by equations (12^ and are not so sensitive to the assumption for 
the viewing angle dependence of r. Since our current procedure looks for the difference with 
respect to the best-fit ringless model, a part of the true signals of the ring is necessarily lost 
by being misidentified as a contribution of the planet. This is why the residual signatures 
of rings shown in Figures E] and [7] are fairly robust against the assumption for the viewing 
angle dependence of r. 



5.4. Detectability of rings 

So far we have ignored observational errors of instrumental noises. In this subsection, 
we discuss the detectability of the ring system and clarify under what conditions the presence 
of rings is statistically confirmed and is really established as an undisputed fact. For this 
purpose, we create mock data and measure the residuals of Avg and F. Using 1,000 mock 
realizations, the statistical significance of the residuals is estimated. 

We define 

No 
i=l 

with the subscript i stands for the sample point. Here, the quantity O denotes either Avg or 
F; aAvs and ap are the statistical errors of radial velocity anomaly and relative flux ratio, 
respectively. We specifically assign the errors of a^vs = 1 ni s~^ and o-p = 10~^ as currently 



Res[C>; 



(24) 
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achieved precision levels, and also of (7^^^ = 0.1m s~^ and ap = 10~^ as future precision. 
The mock data are then created adopting the central values of the fiducial model in Table [1] 
with Gaussian random errors of (T/^^^ and ap. 

In Figure [H we estimate from 1,000 mock realizations and quantify the scatter of 
the functions by varying the number of data samples and the ring parameters of mock 
data. Basically, we choose the central transiting time as the origin of the time and sample 
the data at regular interval. For the cases with future precision (bottom panels of Fig. 9), 
however, we shift the origin of the binning randomly within one bin size with respect to the 
central transiting time. We then explore the effects of bin size by varying the number of data 
points and computing x^- If we do not randomly shift the origin of the binning, shows 
an artificial oscillation due to the sinusoidal nature of the expected signal, and this becomes 
prominent for the cases with higher precision. In order to avoid the artificial coherence of 
the sampled phase, we therefore shift the central position by a size that is randomly chosen 
between to the size of bin. 

Each panel in Figure [H] shows the 95% confidence interval for the distribution of 
values around the mean values {open squares for photometric measurements, crosses for 
spectroscopic measurement). As a reference, we also estimated the scatter of values 
from the mock data without rings, and the resulting 95% confidence interval are plotted 
as solid lines, with middle lines being mean values. Hence, significant deviation from the 
three reference lines disfavors a ringless planet, suggesting the presence of ring as a plausible 
solution. In each panel, the mock data were created adopting the fiducial values listed in 
Table [H except for the obliquity of ring parameters, 6. We specifically examine the three 
cases: 9 = 45° (upper-left), 30° (upper-right), and 60° (lower). 

Figure M indicates that rings with a large obliquity are easily detected. The increased 
number of both the photometric and spectroscopic samples also improves the detectability, 
although the photometric data tend to have a better sensitivity compared to the spectro- 
scopic data. This is partly because characterization of the spectroscopic signal requires two 
additional parameters, i.e., V sinlg and A, in addition to the model parameters of photomet- 
ric signals. Thus, values of the spectroscopic measurement are expected to be generally 
small. Nevertheless, spectroscopic detection/confirmation of rings is valuable and combined 
data analysis will significantly improve the reliability of the detection. This is especially 
relevant for planets orbiting around a rapidly rotating star, in which case the spectroscopic 
signatures of ring during ingress and egress phase become more prominent. 

In Figure [HI we also examine the cases for the close-in planet (HD 209458b) assuming 
a tidally- locked system {left panels) and Saturn system {right panels). The close-in planet 
is likely to be a tidally-locked system and the obliquity of the ring is thus expected to be 



- 20 - 



rather small in such a case. Thus a possible signature of a ring is inevitably weak and it is 
difficult to detect with the current sensitivity (Figure [6]). On the other hand, Saturn has a 
relatively large obliquity {6 = 26.7°) and the photometric signatures of the ring is marginally 
detectable with the HST photometric accuracy of 10~^. The spectroscopic signatures are a 
bit below the detection level of 1 m s~^ partly because the rotation velocity of the Sun is 
small (V^sinJ, = l,880ms-^). 

Nevertheless, the improved precision of future observations will significantly increase 
the detection threshold. Inde ed, the accuracy of spectroscopic mea surement is supposed to 



20081 : iMurphy et al. 1 120071 ) , and Kepler may reach 



achieve ~ a few cms~^ level (ILi et al. 
(Tf ~ 10~^. The lower panels of Figure [9] assume those future precisions. They indicate that if 
possible systematic effects due to the stellar activity are well under control, the spectroscopic 
signatures of rings from the ground-based observation are even detectable. 



6. Conclusions and discussion 

We have explored photometric and spectroscopic detectability of rings around extrasolar 
transiting planets. We have described a general formulation to predict both the photometric 
transiting signature in the stellar light curve and the spectroscopic radial velocity anomaly 
due to the Rossiter effect for a planetary system with a ring. For a face-on ring system, we 
have provided analytic formulae for those effects, which correspond to an improved version of 
our previous results (Paper I). We have shown also several numerical results for more general 
orientations of planetary rings. We found that possible planetary rings around HD 209458b 
and Saturn may be marginally detectable even with the current observational accuracy as 
long as sufficiently sampled data points are available. This is very encouraging, and we 
believe that the detection of rings around extrasolar transiting planets is one of the realistic 
scientific goals of up-coming space missions. 

Let us make a couple of additional comments here. First, our present results may be 
directly applicable to oblate planets (by setting t = oo, Rp = Ri^ = 0, and regarding i?out 
as the surface of the oblate planet, for instance). It would be unlikely that the slow spin 
velocity of tidally-locked close-in planets significantly distorts their sur face, but it would b e 



interesting to put any constraints on the oblateness of the planets (e.g.. ISeager fc Huill2002l ). 



Second, while one woul d expect that the planetary spin axis is aligned with the orbital axis. 



Winn &: Holmad (120051 ) argued an interesting possibility that this is not the case. Thus the 
detection of the ring can test their proposal in a straightforward manner if the planetary 
spin axis is perpendicular to the ring. 
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Fig. 8. — with respect to the the best-fit parameters for a ringless planet model against 
the number of sampled data points (eq.[2l]); the mean value of are plotted in symbols 
with the error bars indicating the 95% confidence interval estimated from 1,000 mock re- 
alizations. Open squares and crosses denote for photometric data and for spectroscopic 
data, respectively The mock data were created adopting the fiducial values listed in Table [H 
except for the obliquities of the ring axis; 6 = 30° (upper- left), 45° (upper-right), and 9 = 60° 
(lower). The spectroscopic and photometric accuracies of (1 m s^^, 10~^) are adopted. 
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Fig. 9. — Same as Figure [HI but for HD 209458 and Saturn systems. Left panels show the 
results for a tidally-locked close-in planet adopting the orbital parameters of HD 209458 
system, while right panels are for Saturn system. The model parameters of these examples 
are described in caption of Figure [61 The upper and lower panels adopt the spectroscopic 
and photometric accuracies of (1 m s~^, 10~^) and (0.1 m s~^, 10~^), respectively. 
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We have to admit, however, that the existence of a Saturnian planet with a ring in extra- 
solar systems is very speculative, in particular for close-in planets where any interpolation 
from our Solar system may not be justified. Thus we share the same feeli ng that Nagaok a 



expressed more than a hundred years ago even in the historical paper (iNagaoka I Il904l ): 
''There are various problems which will possibly be capable of being attacked on the hypothesis 
of a Saturnian system ". Nevertheless we cannot resist the temptation to conclude our paper 
by quoting the last sentence in his paper: ''The rough calculation and rather unpolished 
exposition of various phenomena above sketched may serve as a hint to a more complete 
solution of atomic structure. " 

Finally it would be interesting to note the analogy between the extrasolar planetary 
system and the atomic model; radial velocity modulation of stars due to planets, the Rossiter 
effect of transiting planets, and the signature of rings correspond to the atomic orbital angular 
momentum, the spin of proton, and the spin of electrons. We should not forget that the 
wild idea of Nagaoka led to quantum physics, which turned out to provide the physical basis 
of the astronomical spectroscopic observations via the orbital and spin angular momenta 
of various atomic systems, and therefore he indeed opened a way to the extrasolar planet 
search. 
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A. Approximation of Integrals 

Here, we present the approximate expressions for the integrals Wi {i = 1-4) that appear 
in the formulae fllOj) and fllip . These quantities are all expressed in terms of the planet 
position p = {Xp + Z^y^"^ together with the small parameter 7ps = Rp/Rg. Before writing 
down the approximate formulae, we give precise definitions of these integrals (see Sec. 5 of 
Paper I): 

= [[ dxdz^l - (x2 + z^)/Rl (Al) 

TT-rtp J Js 

W2 = I [ dxdz xVl - (x2 + z^)/Rl (A2) 

XpTCR^ JJs 
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= 11 dxdzVl -x^-P, (A3) 



W4 = / / dxdS iVl-x^-P, (A4) 



where the variables {x, z) means the coordinate system defined in Figure [TJ The variables 
(x, z) are related to the coordinate (x, z) through 



X \ 1 / \ I X 



Note that the integrals Wi and W2 are carried out over the entire planetary disk during 
the complete transit phase, while the regions of the integrals W3 and W4 are restricted to 
the overlapped region between planetary and stellar disks during the ingress and the egress 
phases. 

The approximate expressions for the integrals Wi and W2 are derived by expanding the 
integrand in powers of 7. The resultant expressions become 

Consider next the integrals and W4. They correspond to the ingress and egress 
phases where p ~ 1, and the approximation 1 — 3> 7, which is used in the perturba- 
tive expansion for Wi and W2, breaks down. Thus we attempt to find better approximate 
expressions as follows. Let us integrate first and rewrite W3 and W4 as 

Ws = ^(1 -xo)'(2 + a;o) + / dx g{x;r]p,-f), (A7) 
W4 = -{l-xlf+ di ig{i;r]p,-f). (A8) 



where the function g is given by 

9{x;Vp,V = (l-a;)sm <^ ^— V 



+ ^2{1 + np){xo - - (x - 1 - r,pf}, (A9) 

where 

C . (All) 

2{1 + r]p) 



- 25 - 



Equations (IA7p and flA8|l are still too complicated to integrate analytically. Thus we ap- 



proximate the integrand of equation flA7p as 



9{x- 7) = C, six,, 7) y (,^_,^)(,^_,^_^ + ^) ^ (A12) 
where 

Xc = Xo + siC-l). (A13) 

In the above equations, we determine the two constants s and C3 as follows. From the 
continuity of the first and second expressions of equation (fT2|) . C3 is written in terms of s: 

^(1 -2s7;-7,7) 

The value of s is obtained empirically obtained to yield the best-fit to the numerical results. 
Here we adopt s = 0.4 which gives a slightly better fit than s = 0.5 that was used in Paper 
I. Finally we obtain 



Ws ^ ^(1 - Xo)^(2 + xo) + f 7(7 - W,il - 7). (A15) 

6 2 - 2s7; -7,7) 



We repeat a similar procedure for equation (]A8|) . In this case, the constant C3 is replaced 

by 

C.= ;/f^'' W,il-,), (A16) 

g{l -2s7;-7,7) 

and we adopt s = 0.4 as above. Then an approximate analytic expression for W4 is given as 

^ ^(1 - xir + ^7(7 - 0(2X0 + C - 7) W,{1 - 7). (A17) 

8 4 5f(l - 2s7; -7,7) 
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Table 1. List of parameters 



Variables 


Fiducial value 


Meaning 




Internal parameters of star 




1.146i?o 


Stellar radius 


Ui 


0.2925 


Limb darkening parameter 


U2 


0.3475 


Limb darkening parameter 




Rs^s sin Is = 4700 m s ^ Spin rotation velocity of star 


Is 




Inclination between stellar spin axis and y-axis 


Internal parameters of planet and ring 


Rp 


lM7Rj 


Planetary radius 


A 


0° 


Angle between z-a.xis 






and normal vector fip on {x, 2;)-plane 




1.5Rp 


Radius of inner edge of ring 


-Rout 


2.0Rp 


Radius of outer edge of ring 


r 


1.0 


Optical depth normal to ring surface 


9 


45° 


Obliquity angle of ring axis to z axis 




0° 


Azimuthal angle of ring axis around z axis 


Orbital parameters 


a 


3AU 


Semi-major axis 


-forbit 


1809.01days 


Orbital period 


^ orbit 


90° 


Inclination angle between 






normal direction of orbital plane and y-axis 


b 




Impact parameter, b = (a/Rs) cosiorwt 
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